Partial restoration of factorization and universality in presence of factorization 
breaking interactions in hadronic hard scattering processes. 

A. Bianconl3 

Dipartimento di Chimica e Fisica per I'Ingegneria e per i Materiali, 
Universita di Brescia, 1-25123 Brescia, Italy, and 
Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, 1-27100 Pavia, Italy 

Recent works have discussed the violation of factorization and universality in hadronic hard 
scattering processes aimed at measurements of T-odd distributions. We use simple arguments to 
show that it is possible to restore an approximate factorization involving T-odd contributions, if 
the factorization breaking interactions present a frequency spectrum dominated by a narrow and 
. . . regular peak whose maximum value corresponds to a respected factorization. 
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I. INTRODUCTION 



Recent works [il i] have seriously re-discussed the ideas of factorization and universaUty in hard processes aimed 
at the detection of T-odd distributions. In the last fifteen years the study of these distributions has often made it 
necessary to reconsider the basic properties of the QCD-improved parton model, in particular factorizationj^ 
'" ^ In quark/parton physics the first T-odd distribution function was probably the Sivers function^ in 1990, followed 

I by the Boer-Mulders-Tangerman function^, 0, [S]. The former was used to explain single spin asymmetries Q, the 
latter to explain unpolarized Drell-Yan azimuthal asymmetries [Sj . Recently, it has been demonstrated [To'| that the 
^ I T-odd mechanism introduced in ^11] and ^12,] produces a Sivers asymmetry when its effect is extrapolated into the small 
i ' transverse momentum region. Also, in high energy nuclear physics a T-odd structure function, the so-called "fifth 
structure function", was introduced and modeled p^. [Tsj to describe normal asymmetries in A{e, e'p) quasi-elastic 
scattering. 

For justifying the same existence of leading twist T-odd functions, the central object is the gauge restoring oper- 
ator entering the definition of parton distribution [3]. In an approach where factorization is assumed from the very 
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beginning, this gauge field is assumed to be negligible at leading twist, but in its absence T-odd distributions are 
forbidden IG] by general invariance principles. A chain of arguments and exa mpies[i3, H d IM IM Hi has shown 
that proper taking into account the gauge factor permits the existence of leading twist T-odd distributions in a QCD 
framework, associated with formal factorization breaking by interactions in the initial or final state. In the same works 
' it has been shown that in some cases it is possible to recover factorization from a practical point of view. Systematic 
, efforts have been undertaken ^Ift, ,23J for rewriting factorization rules for processes where partonic distributions depend 
on moderate values of the transverse momentum and the gauge factor is assumed to play a role. 



. !_( Several models and studies have been published on T-odd distributions 124 

together with phenomcnological parameterizations of the Sivers function [3£ 



5-H Tangerman function[8|, l39|. So, we may say that T-odd distributions are an accepted reality, although we have no 
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37, l38|| and of the Boer-Mulders 



281, |29|, |30|, |31|, 133, iSa I3j 



data yet about their behavior at asymptotic energies. 

This creates the problem of the coexistence of factorization breaking and factorization, i.e. at which extent it is 
possible to reconcile the interpretation of observable effects in terms of factorization breaking interactions with the 
safety of the extrernely simple and useful factorization/universality scheme for classifying phenomenology. Also, it 
raises questions [TtI [T9| on the probabilistic interpretation of the measured distributions. Most of the above quoted 
models assume, more or less explicitly, that the factorization breaking effects are there but small enough not to 
overthrow the underlying parton model picture. The analysis of jlj and the observations by |2;] suggest that this point 
is still to be clarified, especially in hadron-hadron single spin asymmetries. 

Without considering a model in detail, we show here that it is possible to restore an approximate factorization 
scheme involving the presence of nonzero T-odd distribution functions, if factorization breaking interactions satisfy 
certain restrictive conditions. 
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The central object of any distribution function calculation is the imaginary part of an amplitude G{2, 1) that 
describes the formation of a quark/hole pair at a spacetime point 1, and the propagation of the hole up to the point 
2 where the hole is filled by the quark again. A value of the longitudinal fraction x is associated with this quark hole. 
By definition x(l) — x{2). In absence of factorization breaking interactions, this is true also in local sense: the quark 
hole conserves its initial longitudinal fraction along all the path from 1 to 2. 

Factorization breaking interactions spoil this property. In any precise model treatment of rescattering, one includes 
explicitly degrees of freedom that exchange 0{Q) amounts of momentum/energy with the quark hole (see e.g. the 
calculations in ref. [l(|), and consequently the wavefunction describing the quark hole contains plane wave components 
associated to values x' ^ x. A way to quantify the effectiveness of rescattering is to look at the size of the difference 
X — x' . 

Our starting point is the idea that a rescattering mechanism produces a finite but small violation of conservation. 
To formalize this, in section IV we introduce the x— frequency spectrum /e(x — x') of the rescattering operator and 
require a set of properties for it. Qualitatively, we may write since now that we require f^ix — x') to consist of a 
regular and narrow peak at a: = x\ with finite width e << 1. In addition, we require fe{x — x') to respect causality 
(see section IV for details). 

The e parameter is supposed to depend slowly on x and on |A:t| (fcr is the quark transverse momentum), and it 
expresses the magnitude of x— nonconservation in rescattering: \x — x'\ e. In the limit where e — > 0, rescattering 
does not affect the quark longitudinal fraction, and factorization is fully respected. 

That the frequency spectrum is large for x' « a;, and small when x and x' are very different, is suggested by the 
fact that at a certain extent factorization works empirically. The assumption that the peak of f^{x — x') at x' = x is 
regular is possibly not justified for individual perturbative diagrams, but becomes more reasonable when we consider 
resummed sets of diagrams, or approximations where some degrees of freedom have been integrated over, or diagrams 
that are regularized via cutoffs, form factors, and so on. Clearly, the chosen values for model parameters may critically 
influence the size of the relevant region of the frequency spectrum. As above written, here we suppose that this size 
e is small. 

If the rescattering operator frequency spectrum satisfies the requested properties, on the ground of simple arguments 
we show that it is possible to write a transverse momentum dependent quark distribution in the form qo + eqi + 
O(e^), where qo may be calculated in complete absence of factorization breaking terms, and qi is a T-odd correction. 

This restores approximate universality for the T-even term: for values of e up to 0.3 it may be rather difficult to 
detect deviations from the qg behavior. This also agrees with the known complete absence of screening effects in 
hadron- nucleus Drell-Yan (see [i^l for a review of the related experiments). 

On the contrary, the T-odd correction term is potentially process-dependent unless one is able to demonstrate that 
e has universal features. An effort in this direction has been undertaken elsewhere [s^] by the author of the present 
work, but the arguments presented here are independent from the analysis in [s^l and so we do not put constraints 
on e in the following. 

As a last remark, the enlarged-delta approximation 15] used in the following was developed to quickly approximate 
interactions between a pointlike particle and a composite system coexisting at the same time. Although it is probably 
possible to re-examine it, and to extend its validity to the case where the interacting particles exist at different 
times (the case of lepton-hadron SIDIS, and interactions between initial and final hadronic systems in hadron-hadron 
SIDIS), the present analysis is limited to the case of Drell-Yan and hadron-hadron initial state interactions. 



II. SOME GENERAL DEFINITIONS AND NOTATIONS 



We work, as far as possible, in the traditional picture of a distribution function, where a quark/hole pair is formed 
in a proton state |P > at a spacetime point "1" = (0, 0, 0, 0), and the hole propagates up to a point "2" = (0, z^, b), 
where the quark is restored with the same momentum it had in the origin. The imaginary part of this amplitude 
gives, by suitable Dirac projections, all the distribution functions we need. 

The quark/hole longitudinal fraction x is defined by P+ = xk^, where P+ and k^ are the large light-cone momentum 
components of the parent proton and of the quark. The infinite momentum limit P+ oo selects leading twist 
contributions. Not to work with singularities, we introduce the scaled longitudinal variable ^: 

C = P+z~, (1) 

and a distribution function may be written as 

q{x,kT) = P+ [ e-"^+'-e'''^''g{z^,kT)dz^dkT = [ e~"^e'''^''G{^,b)d£.dkT. (2) 
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The range of useful values of ^ remains finite ~ when P+ — > oo. So in this limit the function G(^, b) does not 
need to be singular in the origin to produce a nonzero q{x). 

The coordinate z+ plays no role in the following and is not explicitly reported anymore. 



III. ENLARGED 5 FUNCTION 

We adapt here a technique from refs. and [4ll | to roughly but quickly approximate the resummed effect 

of rescattering. Rescattering obliges us to include nonlocality with respect to x, involving integrals of the kind 
/ dx' f(x, x')h{x') where the kernel f(x, x') becomes 5{x — x') in the limit of respected factorization. In this work we 
consider kernels that allow for a small x— nonlocality: \x — x'j ~ e << 1. 

Given a small nonzero parameter e << 1, we define an "e-enlarged (5--function" 5c{x — x') such as to satisfy: 

(i) x' ^ x within uncertainty |e|, 

in) J^^dx'S,ix-x') = l. 
To define St in the limited a;— range [0 — 1], (ii) is better modified into: 

(ii b) 1 — ae^ < J^^ dx'S^ix — x') < 1, with a 1. This property is valid only for x > e and x < 1 — e, so values of 
X near and 1 will not be considered in this work. 

Such functions are normally used to define the ordinary S{x — x') function in the e — > limit. Several functional 
choices for Se{x — x') are possible, the most useful for us is (for e > 0, else e \e\): 

St{x~x') = - . / j-^ 

TT e + [x ~ X 



= 5,+ {x-x') + 5,_[x~x')= ^ \ ^ \ . (3) 

+ ^ ; ' ^ V ; 2iTx-x' + ie 2it x - x' - ie ^' 

Causality implies the substitution 

6,{x-x') ^St±{x-x'). (4) 

that corresponds, in its Fourier transform, to 

expi-e\^\) 9{±OexpiTeO (5) 

Since all the relations interesting us contain a function 9(^) related with causality, Se actually means 5e+, i.e. one 
pole only is encircled in the complex plane integrations. 

We are especially interested in the calculation of integrals of the form dx' f{x')St+{x — x'). For e << 1, and x 
far from its kinematic limits and 1 (a: > e and 1 — a; > e) such integrals may be assumed to be dominated by the 
i5e+ pole x' ~ X ie. 

More precisely, we map the [0, 1] a;— range onto the [— oo, +cx)] z— range by the transformation 

z = ln{x/l — x) (6) 



1 + 



Since dx = e^/(l + e^J^dz, and for x near the complex point a we have l/(a; — a) « (1 + e^)^/e^ ■ l/(z — A) (where 
A = z{a)) we have 

jor X ~ a. (J) 



X — a z — A 

If the integral f{x')Sf^^(x — x') is dominated by the a;' = a;' — ie singularity we have: 



f f{x')5t+{x-x')dx' ^ ^ r — ^^^J^^2^T'^"' 

Jo 27ri [z' ~ z{x + ie)] 



(8) 



f{x + te) = fix) +ie^ + 0(e2). (9) 
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IV. INSERTION OF FACTORIZATION BREAKING INTERACTIONS 



Let us define 

qo{x) = / dfe-"« < Pimi^^mP >, (10) 



q^ix) ^ / dCe-*^« < PimFeiOi^-'mP > ■ (11) 



Of the previous two equations, the former defines a parton-model leading twist quark distribution, the latter includes 
the effect of a rescattering operator i^e(^), where a parameter e quantifying the strength of F is explicitly reported. 
For e = 0, reduces to unity, and q^ to the parton model definition. 

We assume that is scalar, only to simplify the following relations (else, a double sum is present, see the end of 
this section). We introduce the frequency spectrum f^{x — x'): 

F,{0 = I dye-^y^fM, -1 < y < 1. (12) 



For ft{y) we assume the following properties: 

1) Although /e(?/) is not a delta function, it consists of a peak at y = 0. 

2) This peak is regular on the real axis. 

3) The width e of the peak is small: e << 1. 

4) /c(j/) is negligible for \y\ >> e. 

5) Fore ^ 0, Uy) ^ 5{y). 

6) The shape of this frequency spectrum is causality-modified, so that F{£) contains the factor 6'(^). In other words, 
fe{y) is not an ordinary function but a distribution, associated with given integration rules. 

The previous hypotheses allow one to approximate /e [y) the following way: 

Uy) « 5,+ {y). (13) 

Now we may write 

q,{x) = I d^e-"^ < Pm) I dye-'y^S,+ iy)^+iO)\P > = (14) 



dy6,+ {y) J d^e'^'-'-'+y^^ < P|V^(OV'+(0)|P > = (15) 
= dyS,+{y)qo{x + y) = dyS^+ix - y)qo{y). (16) 



Applying the approximation eqlHlwe get 

q,{x) « qoix + ie) « qo{x) + ze^ + 0{e^). (17) 

ax 

If one wants to be more precise and consider the non-scalar operator nature of F^ , one may introduce explicitly the 
matrix sum < P\ ^ ....\P >, and apply the previous considerations to each of the terms of the sum. 



V. DETECTION OF THE FINITE IMAGINARY PART OF THE AMPLITUDE POLE AS A T-ODD 

DISTRIBUTION 

For making the above imaginary shift observable in the form of a T-odd fcy— dependent distribution function, 
factorization breaking interactions must also produce nonzero interference between orbital states differing by one 
unit. 

So, we apply the considerations of the previous subsection to the impact parameter dependent distribution q{x, b). 
As a reference example, we consider the Sivers distribution. A Sivers contribution in agreement with the Trento 
definition (4^ may be introduced in the unpolarized quark distribution of a transversely polarized proton as 
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k 

q{x,k^,0) = quix,k'^) + -^qs{x,k'^), (18) 

where we have assumed that the proton is fully polarized along the y axis. 

We may Fourier-transform this equation with respect to fc, for ky = 0, and write the result in the form 

q{xj) = J dV-''-{Quix,b^) + ib,Q^Qs{x,b^)) = J dV''''-{Quix,b'') + ih,e" (b')Q'!^{x,b')) (19) 

where Qu,s is the Fourier transform of qjj^s with respect to exp{ih ■ fc), and all Q-functions are real and 6— even. 

If we apply the above eq ll7l to the impact parameter dependent distribution, and explicitly separate the 62;— even 
and bx—odd parts of e: 

e = e'ib") + fexe"(62), (20) 

Q,(x,6) « Qo(a; + «e,6') « Qo(x,6) + , /Qo(^'^') (21) 

dx 

= Qo(x,b ) + ie(b ) h ib^e [b ) , (22) 

dx dx 

it is evident that a part of the frequency spectrum of the interaction operator contributes to the Sivers term in eq llSI 

VI. CONCLUSIONS 

In this work, we have considered the action of a set of factorization breaking interactions, whose frequency spectrum 
is regular and peaked near the condition of conserved x. We required the width of this peak to be small: Sx ~ e << 
1. 

We have shown that for such interactions it is possible to retain the factorization scheme, and write the result of 
the calculation of a distribution function as a series expansion with respect to e. 

The zeroth-order term is the T-even and universal distribution function q{x) corresponding to the complete absence 
of initial state interactions. T-even corrections to q{x) appear at second order in e. 

The first order correction is a T-odd distribution function, detectable in angular even-odd interference terms. It 
cannot be stated to be universal, at this level of analysis. 
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